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We present the theory of the leading edge gap in the nor- 
mal state of underdoped high-T c materials. The consideration 
is based on a magnetic scenario for cuprates. We show that 
as doping decreases, the increasing interaction with param- 
agnons gives rise to a near destruction of the Fermi liquid and 
this in turn yields precursors to d— wave pairing. We argue 
that the leading edge gap at ~ 30meV and a broad maximum 
in the spectral function at ~ 150meV are byproducts of the 
same physical phenomenon. 

One of the most intriguing experimental facts about 
underdoped cuprates is that they display superconduct- 
ing properties already at temperatures which can be few 
times larger than T c . This phenomenon, which has been 
observed in the NMR, transport and optical measure- 
ments 0|, is most directly seen in photoemission exper- 
iments on 2212Bi compounds: the leading edge of the 
photoemission curve remains at a finite distance from 
zero energy well above the actual T c and displays an 
angular dependence, similar to that of a true d— wave 
superconducting gap J|,||]. This pseudogap behavior is 
however rather peculiar as the leading edge gap (LEG) is 
not accompanied by the quasiparticle peak. Instead, the 
spectral function is rather flat above the LEG and only 
displays a broad maximum at a frequency ~ 150meV 
which is 5 times larger than the gap. This last fre- 
quency is comparable to the spin exchange integral J, 
and this caused the speculations that the high-frequency 
maximum can be due to the precursors to antiferromag- 
netism A challenging observation for this conjec- 

ture is that the LEG and the broad maximum at higher 
frequencies seem to emerge at the same doping concen- 
tration |g] and therefore are likely to be byproducts of 
the same physical phenomenon. 

In this paper, we show that both, the LEG and the 
broad maximum at higher frequencies can simultaneously 
be explained in the magnetic scenario for cuprates. This 
scenario implies that the low-energy physics of cuprates 
is described by the spin-fermion model in which itiner- 
ant fermions interact with their collective spin degrees of 
freedom by 

U s-f = 9j2 c K^,(3Ck+ q ,pS- q . (1) 

Here g is the coupling constant which is assumed to in- 
crease as the system approaches half-filling, and Ui are 
the Pauli matrices. This model can be obtained from the 
underlying Hubbard-type model by integrating out high- 
energy fermions and performing the RPA summation in 
the particle-hole channel. 



The propagator for low-energy fermions is assumed to 
have a Fermi-liquid form Go{k, to m ) — Zo/{i^ — £k) where 
for e/j we use a tight binding form = — 2t(cosfc x + 
cos k y ) — At' cos k x cos k y — fi consistent with the photoe- 
mission experiments on the overdoped cuprates Q . The 
collective variables S q are characterized by their bare spin 
susceptibility xo(q, w) = x/{l + (q£) 2 - u 2 /A 2 ), where £ 
is the spin correlation length, A = t) s £ _1 where v s is the 
spin- wave velocity, and q = Q — q where Q is either equal 
or very close to the antiferromagnetic momentum (tt, 7r). 
Observe that we did not introduce the damping term 
oc iu> into xo- We argue that for T <C J, the key source 
of spin damping is the decay of a spin fluctuation into a 
particle-hole pair. In this situation, spin damping is not 
an input parameter in the theory, but rather should be 
obtained self-consistently within the spin-fermion model. 

The quasiparticle residue Zq is generally a function 
of T, and it becomes a constant only below some T* [Q, 
when quantum fluctuations start to dominate over classi- 
cal fluctuations. In this paper, we assume for simplicity 
that classical fluctuations can be completely neglected 
below T* and set Zq(T) = Zq{T*). In practical terms, 
this implies that we in fact will be computing the LEG 
right above T c where it is maximal and will not discuss 
how this gap is destroyed by thermal fluctuations. For 
the rest of the paper we absorb both Zq(T*) and x m t° 
the coupling constant: gZo(T*)i/x ~^ 9- 

We now proceed with the calculations. Our strategy 
is the following: we first demonstrate that when the cou- 
pling exceeds some typical value go, the self-energy cor- 
rections to the fermionic propagator nearly completely 
destroy the Fermi liquid in the vicinity of (0, ir) and re- 
lated points. Then we use the renormalized form of G 
to compute the pairing susceptibility in the d x 2_ y 2 chan- 
nel. We show that this susceptibility is attractive, and 
for g > go yields a d— wave LEG above T c . Finally, we 
show how the leading edge gap transforms into a true 
superconducting gap below T c . 

We begin by reviewing the earlier results for the 
fermionic and bosonic self-energies in the spin-fermion 
model [||J|] . The bosonic self-energy gives rise to a damp- 
ing term in the full spin susceptibility: X (g, — 
Xo 1 ^^) + i X~ 1 ujujsf where uj sf « (3/16)v^~ x (g /g) 2 , 
9q = 47ru£/3 and v is the Fermi velocity at the points 
where = ek+Q — t 1 - For the which we are using, 
these points (hot spots) are located near (0, tt) and sym- 
metry related points. 

The fermionic self-energy in the spin-fermion model is 
highly nontrivial even for finite £ due to a hidden singu- 
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larity at ui — > which needs to be regularized, and has 
the form 
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where the two scaling functions have the following lim- 
iting behavior: $i(0) = 1, $(x > 1) — x~ 1/2 , $ 2 (0) = 
1, &2(y > 1) = 41ny/(7ry). Apparently, E(fc,w) is large 
for g ^> go- However, substituting u> s f into the scaling 
functions, we find that for g ^> g , y ^> 1 and hence 
$2(2/) *C 1. In this situation, E(fc,0) depends on g 
only logarithmically and in fact saturates if we impose 
an upper cutoff in the spin susceptibility at u> max ~ 2 J. 
We computed $2(2/) beyond logarithmical accuracy and 
found that (g/go) 2 ®2(y) is always smaller than 1 (it sat- 
urates at about 0.4 for g/go — > 00). In this situation, the 
Fermi surface evolution wouldn't start, and one preserves 
a large, Luttinger-type Fermi surface Q. 

On the other hand, the frequency dependent term in £ 
still scales as (g/go) 2 in a region where x < 1, i.e., where 
e k+Q — 2v 2 (ui/lui) where uj\ = 2w s /£ 2 . In this range, the 
bare ui term becomes overshadowed by the self-energy 
for g > g . Moreover, for e 2 < (9/8)v 2 (lo/uji), the self- 
energy overshadows both u> and e&, and to a good accu- 
racy, the renormalized G acquires a universal, momentum 
independent form 
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where Z = (go/g) 2 . The two conditions on e& and ek+Q 
select a region around a hot spot with the width ~ lu/loi. 
We will see below that the dominant contribution to the 
LEG comes from the frequencies lo ~ u>\. For these fre- 
quencies, Eq. ([}[) is valid over a substantial fraction of 
the Brillouin zone which e.g. includes the (0, 7r) point. 

Eq. (H) is obtained to second order in g but using 
the renormalized form of the spin susceptibility. It turns 
out that higher-order self-energy and vertex correc- 
tions to both, fermionic self-energy and spin damping 
scale in the same way as the momentum dependent term 
in (||), i.e., they depend on g only logarithmically and 
in practice reduce to just constants. We have checked 
that numerically, all third-order corrections are rather 
small and can be safely neglected. For example, for 
(g/go) 2 = 3, the inclusion of the vertex correction into 
the self-energy yields only 4% correction to Eq. (||) [jl0| . 

We now discuss which ratio g/go we expect for cuprates. 
For optimally doped 2212Bi materials, photoemission 
data imply that v « 1.2* ~ OAeV 0. Using £ - 2.5 
inferred from NMR in 214 and 123 materials near T* 
at optimal doping jnj, we obtain uo s f ~ 30(go/g) 2 meV. 
Experimentally, u) s f ~ lOmeV ]ll| which implies that 
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FIG. 1. The photoemission intensity I(u>) = A(u)uf(uj) 
near the Fermi surface for the full solution of the Dyson equa- 
tion with S given by ^ and (g/go)^ = 3 (solid line), and its 
strong coupling version given by (H) (dashed line). We set 
T — 2u>„f. The form of I(ui) is very similar to the measured 
photoemission intensity near (0, n) at optimal doping. 



(g/go) 2 ~ 3. Underdoped materials should have even 
larger ratio of g/go- For (g/go) 2 > 3, we found that 
the full G(k,u) and its strong-coupling version, Eq. (0), 
yield virtually equivalent results for the spectral function 
(see Fig. [|), i.e., for this g/go the spectral weight is al- 
most completely transformed into the incoherent part of 
G, and the quasiparticle peak is hardly visible. 

Notice that though the Green's function in Eq. (||) 
does not give rise to a quasiparticle peak, it has a non- 
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Fermi-liquid form G cx e~ 

\uj\ > uj s f |[2|]. For smaller uj, we have G _1 oc (ui + 
iu>\u)\/(4u) s f)). This is a conventional Fermi-liquid form 
of the fermionic Green's function right at the Fermi sur- 
face. The peculiarity of the present case is that the strong 
self-energy corrections effectively freeze the system at the 
Fermi surface even if actual k deviates from hp and moves 
over the region where Eq. (^) is valid. Away from this 
region, the self-energy corrections get smaller and one 
should recover some renormalized dispersion on a scale 
of t. Notice that this behavior is fully consistent with the 
"flat dispersion" observed near optimal doping . 

Having obtained the form of the quasiparticle Green's 
function, we now consider what happens in the pairing 
channel. The pairing interaction is obtained from (|l|) 
in the same way as in the BCS theory, the only dif- 
ference is that here the intermediate boson is a para- 
magnon rather than a phonon. It has been several times 
in the literature that this interaction is attractive in 
d x 2_ y 2 channel and repulsive in all other channels [fl3"|| . 
The d— wave component of the pairing interaction has 
a form r^,,^ (k, — k,p, — p, ui) — d p T(\uj\) where 
dk = (cos k x cos ky) and r(|^j|) is a decaying function 
of the transferred frequency with the limiting behavior 
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r(0) oc ln£ and T(\uj\) oc lu~ 2 for \lo\ > cj±. Numer- 
ically, r(|w|) is rather flat for \ui\ < lo\ and decreases 
at higher frequencies. To simplify the analysis, we set 
r(|o;|) = r = const for |cj| < lo\ and zero for > lo\. 
The constant is chosen such that the area under T(uj) 
is the same as in the exact expression. This procedure 
yields T = -0.16(g/£) 2 = -3.57 tu sf /Z 2 . Using this ap- 
proximation, we explicitly can sum up RPA series in the 
particle-particle channel and obtain a d x i_ y t pairing sus- 
ceptibility in the form x sc (-fc + q,k,fl) = d\ % sc (g, 0) 
where 



v~fo ft) 3 - r 3rn(g ' 0) 



(4) 



and n(g, f2) = J G(w)G(f2 — w) dpHw is a c?— wave 
polarization operator. The prime to the integration sign 
indicates that the momentum integration goes over the 
region where Eq. ([}]) is valid. In this region, the polar- 
ization operator is independent on q. 

We now compute 11(0, 0) and show that for the fully 
incoherent G from (0), the d— wave pairing interac- 
tion is strongly enhanced, but there is no real insta- 
bility upto T = 0. Indeed, the momentum integra- 
tion in the polarization operator goes over the area 
~ uj/uji. Integrating then over frequency we obtain 
II ~ Z 2 /uj s f f^ 1 (doj/u))(uj/u)i). This integral is clearly 
dominated by lu ~ u>i, and yields II ~ Z 2 /uj s f ~ 
i.e., rn(0, 0) = 0(1) independent on g. Collecting all 
numbers, we obtain 3ITI(0, 0) « —0.7. We see that 
1 + 3rn(0,0) is reduced but still remains positive, i.e., 
fully incoherent G does not give rise to actual supercon- 
ductivity. It does however give rise to d— wave precursors 
as we now show. For this, we construct the pairing self- 
energy using x' sc and obtain 

G-^fc.w) = G- 1 {u) + dl [ x sc (q,tyG(-Lu + ri) (5) 



Here G is the full quasiparticle Green's function, and 
G plays the role of the bare Green's function for the 
Cooper channel. For the 5— functional form of x sc (g, f2), 
Eq. (^|) reduces to a conventional Gorkov's equation for 
the full G. In our case, x sc is enhanced, but it never 
acquires a 5— functional peak. Nevertheless, we can do 
the same trick as with the SDW precursors Q: expand 
G(— ijj + n) = G(— oj) + G' and check whether G' is rele- 
vant. Without G' , (g) has the same form as in the true 
superconducting state. The relative corrections due to G' 
depend the ratio of the typical width of x sc (<7i an( l the 
typical frequency shift obtained by solving the Gorkov's 
equation without G' . If this ratio is small, then the cor- 
rections are also small. Physically, this means that when 
the amount of a shift is larger than the width of a pairing 
susceptibility, the latter can be approximately considered 
as a 5— function at the energy scales comparable to the 
shift. 
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FIG. 2. The photoemission intensity obtained from (^) for 
6 = 3 (dashed line) and b = 0.5 (solid line). The solution for 
b = 3 is valid for T > T c and possesses a LEG and a broad 
maximum at larger frequencies. The solution for b — 0.5 is 
valid for T T c and possesses a quasiparticle peak. As T 
goes below T c , one solution gradually moves into the other. 



We now present the results of computations. Let us 
first neglect G' . Solving (|5|) we then obtain 



G(p) oc 



0J (1 + \J\ - i\u>\ 



-b 2 (l + ^T^J\L 
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where u) — u)/w s f and bk = ZAd(k)/2oj s f where A?(fc) = 
d\ J dqdfl x sc {q,ty- The q— integration in the last for- 
mula again runs over the area where Eq. (Jfy is valid and 
for D, ~ uii which, as we will see, dominates the frequency 
integral, yields 0(1). The integration over fl is not for- 
mally restricted, but the polarization operator decreases 
with increasing such that one progressively looses the 
enhancement in the d— wave channel. We found that the 
polarization operator changes sign at |f2| = u)%. To a 
reasonable accuracy, we can then approximate II(f2) as 
II(fi) = LT(0)(1 - (fi/wi) 2 ). Substituting this result into 
X sc , performing the integration and collecting all num- 
bers, we obtain bk — 0.82£rifc which is a large number 
near (0, 7r) (notice that bk does not depend on g). A 
simple analysis then shows that the pole in G is located 
almost along imaginary frequency axis, at to = —iuj*(k), 
where oJ*(k) — uj s fb 2 — 0.34ri 2 oj\. As a result, the spec- 
tral function which emerges from (|^) does not acquire a 
quasiparticle pole but rather a shift by to ~ to* . For k 
near (0, tt) we then have 



A(uj) OC yful 



■oj*(k) 



uj 2 + (uj*{k)Y 



(7) 



Eqs. (g) and (^) are the key results of the paper. We see 
that already in the normal state the spectral function 
rapidly (as y/uj) increases at low frequencies, reaches half 
a maximum at to rs 0.2a;*, then passes through a maxi- 
mum at lu — to* , and very slowly decreases reaching half a 
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maximum only at u> ~ 5uj* . This behavior has a striking 
resemblance with the LEG behavior observed in photoe- 
mission (see Fig. ||). The position of LEG coincides with 
the half-maximum at low frequencies; the broad maxi- 
mum is located at frequencies which are few times larger. 
This is quite consistent with the data. The magnitude of 
the LEG, iv le « Q.54u) s ^ 2 also has the same order of few 
tens of meV as in the data. 

We now estimate the corrections due to G' . Formally, 
the frequency shift and the width of x sc are both of 
the order of <jj\. However, if the pairing susceptibility 
is strongly enhanced such that 1 + 3rn(0) =i5<1, then 
the width of the pairing susceptibility scales as wi VS. In 
this situation, the relative corrections due to G' scale as 
(5 1 / 2 and are small. In our case, S w 0.3. We computed 
the leading correction to the pairing self-energy due to G' 
and found that near (0, 7r) it accounts for ~ 50% correc- 
tion for u> — uj le , and for only 10% correction for u> = uj* . 
Though corrections are not that small, we expect that 
they somewhat reduce the amplitude of the LEG, but do 
not change substantially the overall shape of A(u>). 

So far we completely neglected the coherent part of G, 
G C oh — Z/{iuj — Ztk)- This piece contributes a conven- 
tional, logarithmical in T term to the polarization oper- 
ator and therefore gives rise to a finite T c . We computed 
T c in a standard manner and found T c ~ vZe~ c ^ where 
c w 1. We see that as the doping decreases, T c actu- 
ally goes down because the correlation length increases. 
Suppose now we are below T c . Then the opening of the 
superconducting gap yields a strong negative feedback 
effect on the spin damping. This gives rise to a rapid 
increase of u) s f compared to the perturbative result, and, 
hence, to a decrease in The latter yields a grad- 
ual shift of the pole in G in (^|) from an imaginary to 
a real axis, which gives rise to a gradual transformation 
of the LEG into the quasiparticle peak (see Fig. ||). 
This is precisely what has been observed in the experi- 
ments dJj|. At T -C T c , bk -C 1, and the typical frequen- 
cies for the pairing problem are much smaller than 0J s f- 
At these frequencies, we have a conventional attractive 
Fermi liquid, which is just frozen at the Fermi surface in 
some A:— range. The pairing then gives rise to a conven- 
tional quasiparticle pole at lu = uj qp = 2bkU s f = ZA d (k) 
(we assume that the total spectral weight in \ sc does 
not change as T goes below T c ). Notice that the posi- 
tion of the pole does not depend on oj s f and hence does 
not change with temperature. We therefore can directly 
compare the locations of oj qp and u le . Substituting the 
numbers, we find co qp — (6/£(T c ))a/ e . For underdoped 
cuprates, £(T C ) ^4 — 5. Then oj qp is slightly larger than 
oj le which fully agrees with the data. 

To summarize, in this paper we have shown that the 
exchange of magnetic fluctuations can account for the 
observed LEG in underdoped cuprates. The LEG and 
the broad maximum of the spectral function at ~ 5 times 
larger frequencies turn out to be byproducts of the same 



physical effect. 

A final point. In the above discussion we assumed 
that the Fermi surface is not modified by interactions. 
In fact, when LEG is formed, the spin damping goes 
down because of a feedback effect from the LEG, and the 
momentum-dependent piece of the self-energy increases. 
Eventually, this increase should trigger the evolution of 
the Fermi surface towards small hole pockets. This evo- 
lution is accompanied by the suppression of the d— wave 
attraction by vertex corrections |^|,[l4| ■ This last effect is 
probably relevant only at very low densities, but it nev- 
ertheless suppresses not only T c but also the LEG before 
the system reaches half-filling. How precisely this hap- 
pens, however, requires further study. 
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